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ABSTRACT 

Detection of Lyman-Break Galaxies (LBGs) at high-redshift can be affected by 
gravitational leasing induced by foreground deflectors not only in galaxy clusters, but 
also in blank fields. We quantify the impact of strong magnification in the samples 
of B, V, i, z &€ Y LBGs (4 < z < 8) observed in the XDF and GOODS/CANDELS 
fields, by investigating the proximity of dropouts to foreground objects. We find that 
^ 6% of bright LBGs {rriHieo < 26) at z ~ 7 have been strongly leased {p > 2) by 
foreground objects. This fraction decreases from ~ 3.5% at z ~ 6 to ~ 1.5% at z ^ 4. 
Since the observed fraction of strongly leased galaxies is a function of the shape of 
the luminosity function (LF), it can be used to derive Schechter parameters, a and 
M*, independently from galaxy number counts. Our magnification bias analysis yields 
Schechter-function parameters in close agreement with those determined from galaxy 
counts albeit with larger uncertainties. Extrapolation of our analysis to z > 8 suggests 
that future surveys with JSWT, WFIRST and EUCLID should find excess LBGs at the 
bright-end, even if there is an intrinsic exponential cutoff of number counts. Finally, 
we highlight how the magnification bias measurement near the detection limit can 
be used as probe of the population of galaxies too faint to be detected. Preliminary 
results using this novel idea suggest that the magnification bias at Muy ~ —18 is not 
as strong as expected if a < —1.7 extends well below the current detection limits in 
the XDF. At face value this implies a flattening of the LF at Muy > —16.5. However, 
selection effects and completeness estimates are difficult to quantify precisely. Thus, 
we do not rule out a steep LF extending to Muy > —15. 
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1 INTRODUCTION 

Surveys of Lyman-Break Galaxies (LBGs) during the epoch 
of reionization aim to make a census of high-redshift galax¬ 
ies and to estimate the available ionizing photon budget for 
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These surveys, however, may provide an increasingly skewed 
view of the early Universe as redshift increases. The ob¬ 
servations are complicated by gravitational lensing ( |Wyithe| 
et al.|20lTl, which affects the observed luminosities and sur¬ 


face density of high-redshift LBGs. Along random lines of 
sight, the probability of significant magnification and mul¬ 
tiple images from gravitational lensing is ~ 0.5% (Barkana 


& Loeb 2000 Gomerford et al. 2002 Wyithe et al. 20111. 


Furthermore, high-redshift luminosity functions have been 
shown to have very steep faint end slopes {a ~ —1.6 at 
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2 ~ 4 to a ~ —2.0 at 2 ~ 8, although the uncertainties are 
large, e.g. Bouwens et al.|2014 Schmidt et al.|[^014 |, which 
results in a bias leading to an enhanced probability of grav¬ 
itational lensing over random lines of sight. This so-called 
magnification bias is further enhanced for flux limits at mag¬ 
nitudes brighter than M* where number counts drop ex¬ 
ponentially. Consequently, bright LBGs become much more 
likely to have been gravitationally lensed than random lines 


of sight (Wyithe et al. 20111. The strongly-lensed fraction 


of LBGs at 2 ~ 7 — 8 brighter than M* is expected to be 

~ 10 %. 

The amplitude of magnihcation bias is a function of 
M* — Miiia and a ( Pei|1995 Wyithe et al.|20fT I, where Mn^ 
is the survey flux limit. Therefore, quantifying the amount 
of magnihcation bias offers a direct probe of the luminosity 
function down to, and below, current survey detection limits 
( Mashian fc Loeb||2013 l. Alternative approaches to quanti¬ 
fying the LF beyond detection limits exist, such as targeting 


use in our analysis. Section describes our method of pre¬ 
scribing a likelihood of lensing to each LBG. Section [^de¬ 
scribes our lensing results, and in Section we assess the 
magnihcation bias and the consequences for the LF beyond 
current survey limits. In Section we assess the observa¬ 
tional ehects of magnihcation bias on the LF and in Sec¬ 
tion we present the strong lensing likelihoods of existing 
2 ~ 9 — 10 LBGs. In Sectionj^we conclude. We refer to the 
HST F435W, F606W, F775W, F850LP, F105W, F125W, 
and F160W bands as 5435, Veoe, * 775 , 2850 , F 105 , ^125 and 
Hieo- Throughout this paper we use Hm = 0.3, Ha = 0.7 
and Ho = 70 km/s/Mpc, and all magnitudes are in the AB 
system ( Oke fc Gunn|19^ l. 


massive galaxy clusters as gravitational lenses (Atek et al. 


2 DATA 

The analysis presented in this paper makes use of the 
wide-area, ultra-deep observations of the XDF/UDF and 


noise characteristics of the background 

(Illingworth et al.|2013 

Windhorst et al.|2011 Grogin et al. 


2011 Koekemoer et al. 

20111. The observations cover the 


Magnihcation bias is expected to signihcantly skew the 
observed bright-end of the luminosity function at very high 
redshifts (see Fig. 3 in Wyithe et al. The effect is in¬ 


triguing in light of the 2 ~ 7 LF, which has been observed 
to both agree with the exponential cutoff in the Schechter 
parametrization (Bouwens et al.||2014|, and also to decline 


less steeply than a Schechter function (Bowler et al.|[2014|. 


4.7 arcmin area of the XDF, which reaches ~ 30 mag at 
5a, the 126 arcmin^ of the GOODS Deep helds, which reach 
~ 28.5 mag at 5 ct, and the 115 arcmin^ of the GOODS 
Wide fields, which reach ~ 27.7 mag at 5a. The catalogues 
were constructed to identify Lyman-Break galaxies from 
2 ~ 4 to 2 ~ 8 using a colour-colour criteria (see details 


Previous studies of high-redshift galaxies have argued that 
gravitational lensing has not significantly affected their lu¬ 
minosity functions (McLure et al.|[2006l, while others have 


made slight corrections in the observed luminosities of LBGs 
due to gravitational lensing (Bowler et al.|[Ml4|. However, 


in Section 3.2.2 of Bouwens et al. 20141. LBGs ‘drop out’ 
in the B 435 , Vooe, * 775 , 2850 and Yios for LBGs at 2 ~ 4, 
2 ~ 5 , 2 ~ 6 , 2~7 and 2 ~ 8 , respectively. The catalogues 
include 5867, 2108, 691, 455 and 155 LBG candidates at 
2 ~ 4 , 2 ~ 5 , 2 ~ 6 , 2~7 and 2 ~ 8 respectively. 

We use the 3D-HST photometric catalogue of the 


even though the sample sizes analyzed in recent studies 
are large, the numbers of bright galaxies, where the mag¬ 
nihcation bias effects will be most apparent, remain small. 
Thus observational verihcation of any changes in slope at 
the bright end remains to be determined. 

Identifying and conhrming individual cases of strong 
gravitational lensing of LBGs at 2 > 4 is made difHcult due 
to sources appearing faint and small. Elongation in the ob¬ 
served LBG due to lensing is difHcult to detect due to their 
small observed size compared with the resolution of the tele¬ 
scope. Secondary images are very difHcult to observe, as they 
will be less magniHed than the primary image, and hence be 
extremely faint ([Barone-Nugent et al.|2013K Secondary im- 


GANDELS area (Skelton et al. 2014 Brammer et al. 


20121 in order to model foreground objects as potential 
gravitational lenses. The 3D-HST survey covers all of 
the GANDELS fields, with spectroscopy compiled from 
the literature. We utilize the spectroscopic redshifts of 
foreground objects when available, and otherwise rely on 


photometric redshifts obtained using EAZY (Brammer et al. 


20081 which are based on deep multiband observations 


(|Erben et aL|2008 

Taniguchi et al.|2007 

Grogin et aL|2011 

Koekemoer et al. 

[2011 

Whitaker et a 

.[ 2011| [McGracken 

et al.|2012 

Bielby et al.j 

2011 [Brammer et al.[[2012[ [Ashby 

et al.[2013 

Barmby et al.|2008[). 


ages will also appear closer to the deHector than the primary 
image, making it more difHcult to disentangle them from the 
deflector galaxy light than the primary image. [Wyithe et al.| 
( |2011| ) calculated the probability of detecting a secondary 
image to be « 10%, where the bright image of a galaxy is 
one magnitude above the survey limit. 

In this paper, we adopt a statistical approach to de¬ 
tect gravitational lensing of high-redshift galaxies using the 
largest samples of LBGs at 4 ^ 2 ^ 8 (Bouwens et al. [20141. 


We calibrate a redshift-dependent Faber-Jackson rela¬ 
tion (Faber & Jackson 19761 in Section [^ using early-type 


galaxy data from Treu et al. (20051, Auger et al. (20091, 


and Newman et al. (20101. The galaxies in these samples 


We assess the likelihood of lensing for each individual LBG 
in homogenous samples at 2 ~ 4, 2 ~ 5, 2 ~ 6 and 2 ~ 7 — 8 
in order to infer the total expected lensed fraction at a range 
of flux limits. In Section[^we describe the data used in our 
analysis. In Section [^ w e derive the Faber-Jackson relation 
(Faber fc Jackson|1976 1 of foreground galaxies that we will 


have published spectroscopic redshifts, velocity dispersions 
and rest-frame 5-band magnitudes. 


3 THE FABER-JACKSON RELATION 

There exists a spatial correlation between bright, high -2 
LBGs and bright foreground objects (see Appendix 0- Spa¬ 
tial correlations between source LBGs and foreground ob¬ 
jects suggests that magnification bias is detectable in current 
surveys. In order to quantify its extent, foreground objects 
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Gravitational Lensing in CANDELS & the XDF 3 


need to be modelled as gravitational lenses. 

The key parameter determining the efficiency of an 
early-type galaxy as a gravitational lens is its velocity disper¬ 


sion (Turner et al. 19841. We estimate the velocity dispersion 


of each galaxy in the CANDELS field from its photometry 
by calibrating a redshift-dependent Faber-Jackson relation 
(Faber & Jackson 1976| FJR). The FJR relates the lumi¬ 
nosity of an object to its velocity dispersion. We include a 
redshift evolution term to account for the evolution of the 
mass-to-light ratio with increasing redshift, so the FJR can 
be expressed as 


Lb = ma2{l + z)^, 


( 1 ) 


where Lb is the R-band luminosity, cr* is the stellar velocity 
dispersion, and 2 is the redshift. The FJR can be expressed 
linearly as 


Mb = ax + hy + c, 


( 2 ) 


where a = -2.57, * = logio 200 kma-i ’ ^ = -2-5^, y = 
log^o(l-I- 2 ), c = -2.51ogi(,(m) and Mb = -2.5 logio(I/s). 
We calibrate the FJR using galaxies with spectroscopic 


redshifts and velocity d ispersions from [Treu et al. (20051, 
I Auger et al.| ( |2009[ ) , and [Newman et al.| ( |2010[ ), v^ch span 
0 < 2 < 1.6. We determine the values of a, b, and c by 
minimising the given by 


X 


-E 


{MBi — axi — hyi — c)^ 
(^Mh. +4t)’ 


( 3 ) 


1=0 ^ Mb 


where are the uncertainties in the data and eint 

is the intrinsic scatter. To avoid degeneracies, we fix the 


slope to be 7 = 3.9, in line with previous studies (Hyde & 
Bernardi|2009a Jonsson et al.|2(il0 l. 


We find best fit parameters of m = 2.3 ± 0.2 x 10® and 
P = 0.7 ± 0.3. The errors on m and /3 are not independent, 
so the uncertainty in the inferred velocity dispersion due to 
their uncertainty is ~ 10 kms“^, and will not significantly 
affect the inferred strongly-lensed fraction. The FJR is 
plotted in the left panel of Figure and the residuals are 
plotted as a function of effective radius, i?o, redshift, 2 , 
and B-band magnitude, Mb are plotted in the centre left, 
centre right and right panels, respectively. The uncertainty 
in the FJR is dominated by the intrinsic scatter, which 
is 46kms“^ in the direction of velocity dispersion. There 
are no systematic biases in the residuals with respect to 
Mb, z or the Re- The scatter in the residuals for galaxies 
at 2 > 0.6 is consistent with galaxies at 2 < 0 . 6 , with no 
evidence of redshift-dependent scatter in our FJR. 

The resultant FJR is consistent with R-band FJRs 
found from weak lensing analyses of type la supernovae 


presented by Jonsson et al. (20101; Kleinheinrich et al. 
( |2004| ); [Hoekstra et al.j ( [2004 1 . As a check of our FJR, we 
compare it with the i*-band FJR (which may be less prone 


to dust-extinction) presented by Bernardi et al. (20031 
for all objects in GOODS. We find very close agreement 
between cr* as inferred from the i-band FJR with our 
R-band FJR. For low-redshift objects, the scatter in the 
residuals between the two methods is 2kms“^. For all 
objects out to 2 = 2 , the scatter in the residuals between 
the two FJRs is 9kms~^. This may be partially due to the 
Bernardi et al. (20031 FJR being calibrated at 2 ~ 0, and 


not taking into account redshift evolution. There are no 


systematic biases in the residuals between these two FJRs 
as a function of 2 , Mb or Re- 

We compare our FJR with estimates of velocity 
dispersions using the stellar mass estimates of galaxies in 
our calibration sample using the relation between stellar 
mass and velocity dispersion presented by [Hyde fc Bernardij 
(2009b I. We find the scatter in the residuals of velocity 
dispersion estimates using this method to be very similar 
(in fact, slightly larger) than those found using our FJR. 
This suggests that using stellar mass information rather 
than Lb will not significantly reduce the scatter in our 
velocity dispersion estimates. 


4 ASSESSING THE STRONG LENSING 
LIKELIHOOD OF LYMAN-BREAK 
GALAXIES 


To quantify the strongly-lensed fraction of LBGs, we model 
every foreground object in the field as a gravitational lens. 
Using photometric information of all foreground objects, we 
ask the following question for each LBG: what is the likeli¬ 
hood of this LBG being gravitationally lensed with magnifi¬ 
cation fi > 2 given its position relative to nearby (in projec¬ 
tion) foreground objects? We disregard deflector-LBG pairs 
with a separation of 6 ^aep > 5'.'0, which is much larger than 
the Einstein Radius of typical deflectors. While the choice of 
maximum separation of 5”0 is somewhat arbitrary, we show 
in Sectionthat 5'.'0 is a reasonable choice. For each fore¬ 
ground object within 5'.ffi of the LBG, we use the following 
process: 

(i) Model the foreground object using a singular isother¬ 
mal sphere (SIS) density profile, 

(ii) Calculate the velocity dispersion that the foreground 
object requires for it to produce an image at the observed 
position of the LBG with a magnification of /r = 2, denoted 
by (J*,req. For an SIS, /r = 2 marks the beginning of the 
strong lensing regime. The required velocity dispersion de¬ 
pends on the LBG-deflector separation, LBG redshift and 
deflector redshift, 

(iii) Calculate the likelihood that the foreground object 
has a velocity dispersion greater than or equal to CT*,req. 
This is the likelihood of strong lensing for that deflector- 
LBG pair, 

(iv) Weight the likelihood of lensing by the inverse of the 
detection completeness at the separation between the LBG 
and the nearby foreground object. 


The final step accounts for reduced sensitivity to faint LBGs 
nearby bright foregrounds. We explain this process further 
in Section (4. II 

To calculate (J*,req we find the Einstein Radius, 0 er, 
required for /r = 2 using the expression for the magnification 
of the image in an observed configuration. 


I ^laep I 

l^sepl — 0ER ’ 


( 4 ) 


where y is the magnification, and 9aep is the observed sepa¬ 
ration between the source image and the deflector. We can 
then find the velocity dispersion corresponding to /r = 2 
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Figure 1. Left: The Faber-Jackson Relation (FJR) that we derive (dashed) projected along the 2 -axis and the galaxies in the three 
samples of |Treu et al.| | |2005[ |, [Auger et ah] l |2009l l; and |Newman et al.|||2010[). We find no systematic biases in the FJR we derive. Centre 
left: The residuals of the velocity dispersions of the galaxies in the |Auger et al.| l |200^ and |Newman et al.| < |20io) | samples as a function 
of effective radius. Centre right: The residuals of the velocity dispersions of the galaxies in the three samples as a function of redshift. 
The scatter in the residuals at 2 > 0.6 agrees with the scatter in the residuals at 2 < 0.6 within the uncertainties. Right: The residuals 
of the velocity dispersions as a function of the R-band magnitude. 


using the expression for the Einstein Radius of an SIS, 

feR = 47r( —(5) 
c Ds 

where a* is the stellar velocity dispersion, Ds is the angular 
diameter distance to the source, and Dls is the angular 
diameter distance from the lens to the source. 

For each LBG-foreground object pair, the likelihood of 
strong lensing of the LBG by the deflector is equal to the 
likelihood that the deflector has a velocity dispersion above 
cr^.req, which is given by 


.if = -erfcl 


u* , 




V^epjR 


)■ 


( 6 ) 


where CT*,inf is the velocity dispersion inferred from photom¬ 
etry (using the FJR), and cfjr is the intrinsic scatter in the 
velocity dispersion of the FJR. 

In the event that there are multiple potential deflec¬ 
tors within 5f0 of the source, we treat them independently 
and calculate the probability that at least one is lensing the 
source by /r ^ 2. For n deflectors, this is 


n 

^=i_n(i-if,). (7) 

We show a subset of the sample consisting of some of 
the highest-likelihood lenses in Figure We describe these 
systems further in Section [5.1| 


4.1 Accounting for Sensitivity Variations 

Faint LBG samples have a reduced completeness compared 
to bright ones. This alone would not affect our inference of 
the lensed fraction, because the completeness would change 
the numerator and the denominator by the same factor at 
fixed magnitude. However, we note that there may be a 
further reduced sensitivity to detecting faint LBGs around 
bright objects, which will affect potentially lensed LBGs dif¬ 
ferently to those isolated in the field. This effect could cause 
our measured strongly lensed fraction to be artificially low. 

We weight all LBGs that appear close in projection to 
bright foreground objects by the inverse of their relative de¬ 
tection probability in order to account for reduced sensitivity 


around foreground objects. To do this, we run completeness 
simulations around all foreground objects which are either, 

(i) Assessed as having a greater than 1% chance of lensing 
a nearby LBG, or, 

(ii) Brighter than = 24 mag and within 2"5 of an 
LBG. 

We run source recovery simulations in order to determine 
completeness as a function of radius around each foreground 
meeting either of the above criteria. The source recovery 
simulations are run for LBGfS at the redshift and of the 
magnitude corresponding to that of the nearby LBG. The 
completeness of a source LBG, Sc, becomes unaffected by 
typical foregrounds at a separation of around l"5 —2'.'0. The 
weight, Wc, we apply to each LBG is defined as the inverse 
of the completeness, Wc = l/sc- We apply a maximum 
weighting of Wc. = 10 to any LBG. 

We find that weighting the lensed fraction in this way 
has a minimal effect on the bright-end of the observed 
lensed fraction. However, the relative completeness of faint 
galaxies near to bright foreground galaxies is, as expected, 
lower than for the brighter LBGs.. 


4.2 Uncertainty Checks 

The typical uncertainty in the photometric redshifts of 
foreground sources derived in the 3D-HST catalogues are 
/S.Z ~ 0.1 — 0.2, so the uncertainty in the magnification 
is dominated by the intrinsic uncertainty in the FJR. 
Furthermore, the uncertainty in Mb due to photometric 
redshift errors (via the distance modulus) is partially self¬ 
regulating as the inferred velocity dispersion is a decaying 
function of redshift, while inferred rest-frame luminosity is 
an increasing function. We find that of the 40 2 ~ 7 LBGs 


^ The artifical sources in the recovery simulations are extended, 
and have sizes typical of LBGs at the appropriate redshift. It 
should be noted that strong lensing may cause the sources to 
appear more extended, which will affect their completeness. This 
effect is expected to be small, but is a slight limitation of this 
analysis. 
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Figure 2. Examples of possibly-lensed LBGs in the four samples. All cutouts are of the J 125 images, are 10”0 X 10”0 and are shown on 
the same contrast scale (except for the top, left cutout, which contains two very bright foreground galaxies). The LBGs are circled in 
red and the deflectors are labelled by their spectroscopic/photometric redshifts. Each LBG is labelled with its L/igo magnitude, and its 
likelihood of being strongly lensed (top left corner). The LBG shown in the bottom, left panel is the brightest LBG in the z ~ 7 sample. 


that have a likelihood of lensing of ^ 10%, the deflector of 
only one has a photometric redshift with less than an 80% 
chance of residing within Az = 0.2. None of the deflectors 
of z ~ 4 — 6 LBGs with a likelihood of lensing of 10% 
have photometric redshifts with less than an 80% chance of 
residing within Az = 0.2. The uncertainty in source redshift 
(Az ~ 0.35) is negligible as the angular diameter distance 
is a relatively flat function at high-redshift. 

There is a possible Eddington bias stemming from 
uncertainties in the photometry and the shape of the LF 
at z ~ 1, which could bias the inference of cr* from the 
B-band luminosity. We find that 'I'(L ± 5L) only varies by 
2 — 5% from <['(1/) for galaxies brighter than ~ M* at z = 1 
in the field. Therefore the number density of bright galaxies 
does not change significantly within the photometric 
uncertainties, and the Eddington bias is negligible. 

We note that a limitation of this analysis is that it 
assumes all deflectors to be singular isothermal spheres. 
Singular isothermal ellipsoids (SIEs) may be a more realistic 


parametrization of potential deflectors (see Keeton 2001 


for ellipsoidal density parametrizations). However, SIEs do 


complicate the calculation signihcantly (Kormann et al. 


1994 Huterer et al. 20051. The median ellipticity for all 


objects in the CANDELS fields is e = 0.21, with 83% 
having an ellipticity of e < 0.4. In the case of low deflector 
ellipticity (e < 0.2), the change in the magnification 

estimate is ~ 10% along both the major (+10%) and 
minor (—10%) axes. For larger ellipticities (e ~ 0.4), the 
magnification estimate becomes « 20% lower for an image 
located along the minor axis, and a factor of two higher for 
images along the major axis. Using an elliptical deflector 
model for the system shown in the bottom, centre-left panel 
of Figure which includes a deflector with large ellipticity 


(e = 0.48), we estimate the magnihcation to be /r ~ 1.6, 
as compared with the SIS estimation of /r ~ 1.8. The LBG 
in the bottom, right panel of Figure which is near a 
deflector with ellipticity e = 0.33, has a magnification of 
/r ~ 1.6 in the SIS model, which becomes fj. ~ 1.9 using an 
ellipsoidal model. 

Similarly, the lensing cross section (the strong lensing 
area in the image plane owing to a deflector) of an SIE is 
the same as the optical depth of a SIS with a higher-order 
term (Kormann et al. 19941. The area of sky covered by the 
Einstein Radius of an SIE is only ~ 5% larger than the area 
of sky covered by an SIS for reasonable ellipticities (e < 0.4). 
Therefore our calculations of the optical depth in Section 
are not significantly affected by the SIS assumption. 
These calculations are consistent with previous studies of 
the effect of ellipticity on the strong lensing optical depth 
and magnification, such as Huterer et al. 120051 who noted 
that aside from image multiplicities, introducing shear and 
ellipticity has surprisingly little effect. Hence, using an SIE 
deflector will not qualitatively change either our strongly 
lensed fraction or magnification bias results. 


5 THE STRONGLY LENSED FRACTION 

The method of prescribing a likelihood of strong lensing 
described in Section was applied to each LBG in the 
samples at z ~ 4, z ~ 5, z ~ 6 and z ~ 7 — 8. The 455 zsso 
dropouts and 155 Yios dropouts were combined to create 
a statistically significant sample with a mean redshift of 
z = 7.2. The number of strongly-lensed LBGs brighter than 
a theoretical survey limit, mum, is assessed for each of the 
samples. For the 1775 and zgso & Ylos samples, we assess the 
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Figure 3. Left: The lensed fraction of background LBGs as a function of flux limit for the 2850 Tios-dropout samples (blue), 
2775 -dropouts (red, offset by m -h 0.1), Veoe-dropouts (green, offset by m -h 0.2) and B 435 -dropouts (yellow, offset by m -h 0.3). The 
observed lensed fraction decreases monotonically with decreasing redshift for flux limits of miim = 26, 27, 28 30. The analytic strong 

lensing optical depths, r, (strongly-lensed fraction of random lines of sight) for each source redshift are plotted as dashed lines in the same 
colours as the observed lensed fractions (see Section Right: The observed magnification bias at each redshift overlaid as a function 
of M* — Miini (see Section |^. The -B 435 sample contains the brightest measurements with respect to M* (1 mag brighter), followed 
by the Vgoe sample, the 2350 & ^105 sample, and the *775 sample. The bias is defined as the ratio of the solid and dashed lines in the 
left panel. We show the bias at each redshift individually in Figure]^ For an LF without strong evolution in a, which is approximately 
observed from 4 < 2 < 7, the bias is not expected to evolve. A roughly constant bias is observed at all values of M* — Mum for the four 
independent LBG samples from 4 < 2 < 7. 


lensed fraction brighter than mum = 26, 27, 28, 29 & 30 mag 
in 7 /i6o- We include mum = 25 mag for the Vgoe sample 
and miim = 24 & 25 for the B 435 sample, as M* appears 
brighter for these samples. The strongly lensed fraction is 
not affected by the differing depth of the CANDELS fields 
and the XDF. 

The cumulative lensed fraction at each of these flux 
limits is the ratio of the expected number of strongly lensed 
LBGs (the sum of all lens likelihoods) brighter than the 
flux limit and the total number of LBGs appearing brighter 
than the flux limit. However, the cumulative lensed fraction 
depends on the total completeness of the combined sample. 
To account for incompleteness in number counts, we use 
the LF of [Bouwens et ^ (20141. The cumulative lensed 
fraction at 2 ~ 4, 2 ~ 5, 2 ~ 6 and 2 ~ 7 is shown in the 
left panel of Figure]^ The right panel of Figure]^ shows the 
observed magnification bias (see Section]^. When inferring 
properties of the luminosity function (Sections |6.1| & |6.2[ ) 
we use the observed lensed fraction in each bin without 
LF-correction so as to not presuppose the nature of the LF. 

The trend to a larger fraction of strongly-lensed galax¬ 
ies for brighter flux limits is reasonably smooth, monotonic 
and observed in each of the four independent samples. 
The amplitude at all flux limits steadily increases from 
2 ~ 4 to 2 ~ 7 (although the error bars are large), which 
is expected as the faint-end slope steepens and the strong 
lensing optical depth increases at higher redshift. The 
excess probability of gravitational lensing of bright galaxies 
is detected at high-significance in each of the samples. The 
lensed fraction of LBGs brighter than rriHigo = 26 is ~ 6% 
at 2 ~ 7 and ~ 3.5% at 2 ~ 6, although the uncertainty is 
large due to the rarity of bright objects at high-redshift. At 


2 ~ 5, the lensed fraction at the same flux limit is ~ 3.5% 
and at 2 ~ 4 the lensed fraction is ~ 1.5%. 

We also assess the lensed fraction at brighter flux limits 
for the 2 ~ 4 and 2 ~ 5 samples. We find that the lensed 
fraction continues to rise, as expected. At 2 ~ 5, ~ 5% of 
LBGs brighter than m = 25 are strongly lensed, and at 
z 4:, ~ 4.5% of LBGs brighter than m = 24 are lensed. 

The errors are calculated using bootstrap resampling. 
The bootstrap sample is drawn from the entire sample 
with replacement N = 10^ times. Each time, each LBG 
is considered either “lensed” or “not lensed” randomly 
according to its likelihood of having been lensed. The lens 
fraction is recalculated for all limiting fluxes. The error bars 
represent the Icr limits of the resultant distributions. 


5.1 Examples of Likely Lensed Systems 

We present an illustrative sample of some llkely-lensed 
candidates in the surveys in Figure Gases from our 
highest -2 sample are emphasised because they are of the 
most interest, and have the most importance to future 
surveys. We note that the three brightest 2 and T-dropouts 
in the entire sample are each deemed to have a likelihood of 
lensing of > 10%. The brightest LBG in the 2 ~ 7 sample 
is shown in the bottom left panel of Figure All cutouts 
are shown at the same contrast scale, except for the 2 ~ 4 
lens candidate (top left), which is in proximity to two very 
bright foreground galaxies, both with Mb ~ —23.5, one of 
which is spectroscopically confirmed at 2 = 0.8. All cutouts 
are 10”0 on each side. In each case, the deflector candidate 
is labelled with its spectroscopic or photometric redshift 
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and the LBG is labelled with its iifieo magnitude. 

The cutouts highlight the difficulty in locating sec¬ 
ondary images in the event the LBG has been strongly 
lensed. A secondary image will appear closer to the fore¬ 
ground galaxy than the primary (circled) image, and is 
likely to also be appear much fainter than the primary 
image. 


5.2 Deflector Properties 

We present the distribution of the image-deflector sepa¬ 
rations, deflector redshifts and deflector B-band absolute 
magnitudes in this section. The number of lensed sources 
is weighted by the likelihood of lensing for each image- 
deflector configuration. 

The top row of Figure shows the distribution of 
lens rest-frame B-band magnitudes for each of the four 
independent LBG samples. The peak of the distribution 
occurs around Mb ~ —22 for each of the samples. 

The middle row of Figure shows the distribution of 
image-deflector separations for each of the LBG samples. 
The normalised cumulative fractions are shown as dashed 
lines. We observe an approximate increase in the peak of 
the separation distribution as redshift increases (from ~ I'.'O 
at 2 ~ 4 to ~ 2'.'0 at 2 ~ 7), consistent with the expectation 
that higher-redshift sources have larger deflection angles. 

The bottom row of Figure shows the distribution of 
deflector redshifts. The normalised cumulative fractions are 
shown as dashed lines. We observe an increase in the peak 
of the deflector redshift distribution from 2 ~ 4 sources, 
where the deflector distribution peaks around 2 ~ 1, to the 
2 ~ 7 sources where the peak occurs around 2 ~ 2. This 
evolution is consistent with the expectation that lenses 
are most likely to be found at around half of the angular 
diameter distance to the source. 


6 MAGNIFICATION BIAS 


The total magnification bias of a flux-limited sample, B{L > 
Liim), is the ratio of the fraction of strongly-lensed galaxies 
and the fraction of strongly-lensed random lines of sight, de- 


hned as the strong lensing optical depth, r (e.g. Wyithe et al. 
|2011[ ). We generate a catalogue of 50,000 random source po¬ 
sitions in the GOODS fields and use the method of assessing 
the lensed fraction presented above in Section to deter¬ 
mine the fraction of the source plane that will be strongly 
lensed. Based on our FJR, we assess the strong lensing op¬ 
tical depth for sources at 2 ~ 4, 2 ~ 5, 2 ~ 6 and 2 ~ 7.2 to 
be T = 0.41%, 0.54%, 0.65%, 0.75%. The values found are 
broadly consistent with theoretical predictions of the strong 
lensing optical depths at these redshifts (Mason et al.|20l5 


Wyithe et al. Mil [Barkana fc Loeb|2000| . Due to the large 

number of foregrounds in the CANDELS fields, the relative 
statistical uncertainty on r is only ~ 4% in all samples, and 
hence negligible in our bias calculations. We And consistent 
values for the optical depth if we apply a reasonable upper 
limit of ~ 350kms“^ on the inferred velocity dispersion of 
foreground galaxies. We find that our method of determin¬ 
ing the optical depth returns values in close agreement with 


those in Mason et al. (2015 I when we adopt their method of 
inferring velocity dispersions using stellar mass estimate^ 
The optical depths are plotted as dashed lines in Figure!^ 

The bias is therefore the observed magnified fraction 
divided by the optical depth (the solid lines divided by the 
dashed lines in Figure]^. The observed total bias for each 
of the samples at a range of flux limits is plotted in the 
right panel of Figure and the top row of Figure The 
bias reaches values of ~ 10 at bright magnitudes and high- 
redshifts, but near the survey flux limit has values lower 
than expected for a LF that remains steep well beyond sur¬ 
vey limits. 

We calculate the observed magnification bias in each 
bin (as opposed to the total magnification bias for all galax¬ 
ies brighter than a flux limit). The results are plotted in the 
bottom row of Figure 

For a LF with weak (or no) redshift-evolution of the a 
parameter, the magnification bias as a function of M* — Miim 
is expected to remain approximately constant with redshift. 
To highlight that this trend exists in the data, we plot the 
observed magnification bias at each redshift on the same 
axes in the right panel of Figure While a evolves from 
~ —1.6 to ~ —2.0 from 2 ~ 4 to 2 ~ 8, the statistical un¬ 
certainties in our measurements are larger than the change 
in bias from this evolution. 

For a given luminosity function, 4'(L), the magnifica¬ 
tion bias can be predicted anal3rtically ( [Turner et al.|[T984[ ) 
at luminosity L, by 


B(B) = 


T(L) 


( 8 ) 


where the 1/p factor accounts for the stretching of dfi with 
magnification, and ^ is the magnification distribution for 
the brighter image in a strongly lensed system, given for an 
SIS by. 


dfi \ 0 


for 2 < p < 00 
for p < 2 


(9) 


We assume 4/(L) to be the Schechter luminosity function. 
The analytic magnification bias for all galaxies in a flux 
limited sample is 


B(L > Blini) — 


I, 


A^max 


dmL) 


( 10 ) 


where the factor of 1/p is no longer included because we 
are now integrating over magnification, and the stretching 
of dp no longer affects the integral. 

Results for our bias estimates given the Schechter 


LF parameters in Bouwens et al. (20141 and theoretical 
curves are plotted in Figure [^ The top row compares the 
theoretical bias of all galaxies in a flux limited sample with 
our measurements of the observed bias for all galaxies in a 
flux limited sample. The bottom row shows the theoretical 
bias of galaxies at a fixed luminosity with our measurements 
of the observed bias in each magnitude bin. Theoretical 
values for bias are calculated using previously derived LF 
parameters a and M* (Bouwens et al. 20141 and a range 


^ We use stellar mass estimates of foreground galaxies in the 
GOODS/CANDELS fields from the 3D-HST catalogue 
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z=3.8 



z=4.9 


z = 5.9 


z = 7.2 












Figure 4. The lensing-likelihood weighted distributions (solid) and cumulative distributions (dashed) of deflector properties. These 
distributions illustrate the diversity and evolution of the deflector population in the four samples analyzed. Top row: The distribution 
of S-band absolute magnitudes of the deflectors for the four LBG samples. Middle row: The distribution of image-deflector separations 
for the four LBG samples. Bottom row: The distribution of redshifts of the deflectors for the four LBG samples. 


of values at which the luminosity function deviates from a 
steep faint-end slope. We And close agreement between the 
observed shape and amplitude of the magnification bias and 
the theoretical function in each of the independent samples. 

It is worth noting that the inferred magnification 
bias is not sensitive to the parameters of the Faber- 
Jackson relation (Section [^, because the use of the FJR 
to determine the efficiency of observed galaxies affects 
both the numerator (fraction of strongly-lensed LBGs) and 
the denominator (the strong lensing optical depth) similarly. 


6.1 The Faint-end Slope Beyond Current Flux 
Limits 


Magnification bias results from magnification of intrinsically 
faint sources below an observed flux limit into an observed 
sample, hence quantifying the degree of magnification bias 
offers an opportunity to investigate the behaviour of the LF 
beyond current survey limits. 

To illustrate, we begin with a toy model in which there 
is a minimum luminosity for galaxies of Lmin, below which 
there are no galaxies, and a power-law slope of a = —2.0 for 
L > Lmin- In this toy model, a sharp cutoff in the LF at a 
value of Lmin yields a bias of. 


Bi^L > Llim) 


3 j It ■ 1 -^lim ^ 2Z/rnin 

1 for Lmin < Llim < 2Lmin 

( 11 ) 


This implies the total bias of a flux-limited sample reaches 
unity approximately 1 mag brighter than Lmin- We observe 
a hint of the possibility of this occurring in the four samples 
presented in this paper, as seen in the top row of Figure]^ 
The bottom row of Figure also highlights this behaviour 
in our samples. 

Rather than a sharp cutoff in the LF, we consider a more 
realistic model with a LF that flattens (02 = — 1) after some 
luminosity, Ltum. We attempt to constrain Ltum by hnding 
a LF that will reproduce the observed magnification bias 
in the bottom row of Figure Using a broken Schechter 
function of the form. 


^{L)dL = 



for L ^ Ltum 

for L < Lturn 
( 12 ) 


with 012 = —1 (i.e. a flat LF beyond Ltum), we ht the bias 
calculated from equationto the data with Ltum, cti, and 
M* as free parameters. We include priors on the values 
of ai and M*, from the LFs of Bouwens et al. (20141. 
Figure shows the constrains found for the minimum 
luminosity from our analysis fitting to two subsets of our 
measurements. The dashed black line in Figureshows the 
probability distribution function (PDF) when fitting Mtum 
to the observed bias of only LBGs brighter than m = 29, 
and the solid line shows the PDF when fitting to LBGs 
brighter than m = 30 (approximately the 5cr limit in the 
XDF). The PDFs are normalised such that the probability 
of Mturu < —12 is unity. We find preferred values of Mtum 
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z=3.8 z=4.9 z=5.9 z = 1.2 





rrijj rrijj m^j rrijj 


Figure 5. Top row: The observed total magnification bias of all galaxies brighter than a flux limit M* — (solid) compared with 

the theoretical magnification bias for the luminosity functions (described by equation |10[ l from [Bouwens et ah] l |2014[ | with a range 
of magnitudes at which the LF flattens (a ~ —1), denoted by Mtum- The cumulative lensed fraction is corrected for incompleteness 
according to the LF of [Bouwens et ar] ( |2014| |. Bottom row: The observed magnification bias in each bin plotted at the mean luminosity 
of the bin. At faint magnitudes in each sample of LBGs, the bias falls below the value expected from a faint-end slope continuing well- 
beyond the survey limit, indicating a possible deviation from a steep faint-end slope of the LF, although they agree with theory within 
their error bars. For a LF with a steep faint-end slope continuing well beyond the flux limit, the bias flattens to a value of B ~ 2 — 3 
(depending on a). The measurements of bias at a fixed luminosities do not need to be corrected for incompleteness. 
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Figure 6. The inferred value of Mtum using the magnification bias measurements of all galaxies brighter than m < 30 (solid) and all 
galaxies brighter than m < 29 (dashed). The bias measurements which we fit to are shown in the bottom row of Figure]^ In black we 
plot the likelihoods with a prior on a and M* from [Bouwens et al.j | [2014| | (see Table 4 therein for values). The red curves show the 
estimated likelihoods including an additional requirement that the minimum magnitude is fainter than the magnitude to which current 
observations confirm a steep faint-end slope. We find approximately consistent preferred values of Mtum in each of the four samples, 
but with varying amplitudes. For LBGs brighter than m = 30 (approximately the 5cr XDF limit, solid line), we find a preferred value of 
Mtum around the observational limits in each sample. However, a value of Mtum > —16 is not excluded. For only LBGs brighter than 
m = 29 (dashed line), we find no constraint on Mtum in any of the samples. This is expected, because to constrain Mtum we need to 
consider galaxies within ~ 1 — 2 mag of Mtum- The curves are normalised such that the probability of Mtum < —12 is unity. 
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to peak around the current observational limits in each 
of the independent samples from 4 < 2 < 7. The sample 
at 2 ~ 6 peaks at a magnitude brighter than flux limits, 
which can be ruled out observationally. We also calculate 
the likelihoods with an additional prior enforcing Mtum 
to occur below the magnitude that the steep faint-end 
slope has been observed to extend to. These are plotted 
in red in Figure for the same flux limits as above. It is 
important to note that the inference of Mtum occurring 
close to current flux limits is marginal and does not rule out 
a faint-end slope extending well-beyond current flux limits, 
or a flattening for a few magnitudes followed by an upturn. 

We find that the constraint disappears when Mtum is 
fit to only brighter (m < 29) galaxies. We extend this test 
by recalculating our results by omitting XDF LBGs entirely 
from the analysis to investigate whether the observed 
magnification bias will always approach unity near the flux 
limit due to selection effects. When we perform this test, we 
find the magnification bias of galaxies in the GOODS-North 
and GOODS-South is completely consistent with that of 
the full sample, rather than approaching unity near the flux 
limit. 

It is important to note that the magnification bias 
is only observed to drop below its expected value close 
to the current flux limits where selection effects become 
significant. While we have taken care to account for the 
decreased sensitivity to very faint sources, there still exists 
the possibility that we have missed a significant fraction 
of gravitationally lensed LBGs at very faint magnitudes. 
Furthermore, if the interloper fraction at very faint fluxes 
is high, the lensed fraction will be underestimated, causing 
a spurious inference of Mtum- 

The possibility of a flattening of the LF at 
Mtum ~ —16.5 at 2 ~ 7 is consistent with observa¬ 
tions of LBGs down to Muv ~ —15.5 of magnified 2 ~ 7 
LBGs using Frontier Fields cluster Abell 2744 ( |Atek et al. 
2014 Ishigaki et al. 2014[). F igure shows the [Bouwens 
et al. (20141 and Atek et al. ( 2014[ ) data with the best fit 
broken Schecter function (02 = — 1). We find that a broken 
Schechter function represents the data very well, and offers 
an independent constraint on Mtum- The data favours a 
value of Mtum ~ —16.5, consistent with Mtum inferred from 
our magnification bias results. However, as is the case with 
our magnification bias analysis, this inference is based on a 
single data point. 

Further to this study at 2 ~ 7, there exist measure¬ 
ments of the UV luminosity function of 2 ~ 2 LBGs down 
to Muv ~ —13 ( |Alavi et al.|2014| ). While a single Schechter 
function favours a steep slope extending to Muv ~ —13 
at 2 ~ 2, the shape of the LF at the faint end may be 
more complicated than this simple parametrisation. In 
fact, the data also seems to suggest a flattening of the 
population density between —17 < Muv ^ ~15, before a 
steeper upturn from —15 < Muv ~13 (see |Alavi et al.| 
2014 Figure 7). This more complicated LF would produce 
a reduced magnification bias {B{L > I/um) ~ 1) of 2 ~ 2 
LBGs near a flux limit of Muv ~ —17.5, which is on the 
edge of current survey limits in blank fields ( Sawicki||2012 


Oesch et al.|2010 Hathi et al.|20l6 Reddy fc Steidel|2009 1. 
Additionally, a flattening, and potentially a rise in the LF 
at Muv ^ “15 would not be inconsistent with the inference 


from GRB host galaxies studies (Trenti et al. 2012 2013 



0.5 


0.4 



0.1 


0.0 


Figure 7. The 2 ~ 7 LF data from [Bouwens et al.| | |2014l l (blue 
circles) and jAtek et al.||2014| | (black squares) with the best fit bro¬ 
ken Schechter LF (dashed), and single Schechter LF (dotted). The 
solid line shows the likelihood of Mtum given the two datasets. We 
find a favoured flattening magnitude at Muv ~ —16.5, consistent 
with our magnification bias measurements. 


Tanvir et al. 20121. In fact, these studies only constrain 


the presence of an abundant population of galaxies below 
the XDF detection limit, but not the shape of the galaxy 
LF, which they assume to be Schechter-like. Interestingly, 
theoretical models that are based on a double population 
of faint galaxies have been proposed in the context of 
hydrogen reionization (e.g., see Alvarez et al.|[2012|. 


6.2 Deriving Schechter Parameters from Lensing 

A very interesting application of this analysis is that 
Schechter function parameters a and M* can be derived 
directly from the magnification bias. This method is 
completely independent of the standard procedure using 
number counts of galaxies, and therefore could be combined 
to produce improved constraints. 

The magnification bias at a fixed luminosity can be 
predicted using equation and is a function of a and M*. 
By fitting the predicted bias of LBGs at a fixed flux to 
our measured bias in each flux bin (which is shown in the 
bottom row of Figure and is not the LF-corrected cumu¬ 
lative fraction), we can constrain a and M-!^ This does not 
rely on any prior knowledge of the LF. Because a and M* 
are much more sensitive to the bias of bright galaxies than 
that of faint galaxies, and we see a possible deviation from 
a single Schechter function at faint magnitudes, we exclude 
the two faintest bins ^9 < m < 30, and 28 < m < 29) 
from the fil|^ Figure shows the constraints on the LF 
from the observed magnification bias alone along with the 
constraints from number counts of the same samples. 


^ We fit the theoretical bias at the mean magnitude of the LBGs 
in each magnitude bin to the observed bias in that bin with a 
and Mi, as free parameters 

In fact, fitting the data including the two faintest bins with an 
extra free parameter, Mturm a^nd marginalising over this param¬ 
eter gives the same result 
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z = 3.8 


z = 4.9 


z = 5.9 


z = 1.2 














Figure 8. Measurement of the Schechter parameters, a and M* using only the observed magnification bias as a function of M* — 
(black). The contours from [Bouwens et al.| | [2bl4[ | (red) and [Finkelstein et al.| ( [2014| | (yellow) are shown for comparison. We find close 
agreement between the two methods at 2 ~ 4 and 2 ~ 5, while the constraints at 2 ~ 6 and 2 ~ 7 from lensing are weaker due to the 
larger error bars on the magnification bias measurements. 


The measurements of the bias at 2 ~ 4 and 2 ~ 5 do 
an excellent job of constraining the luminosity function. 
At higher redshift, as the samples become smaller and 
the random errors grow, we cannot constrain the LF as 
effectively. However we find that our observations are 
consistent with the UV luminosity functions presented 


by Bouwens et al. (20141. This also provides an internal 


consistency check of our analysis. 


6.3 Contaminant Discussion 

We check for bias arising from the selection of LBGs. There 
may be an enhancement and reddening of LBG candidates 
observed around bright, red foreground galaxies due to 
photometric scatter, causing an increased fraction of inter¬ 
lopers around such foreground objects and providing a false 
lensing signal among bright candidates. To determine if this 
may affect our results, we check if an enhanced interloper 
fraction around bright foregrounds is found in lower-redshift 
LBG samples for which there exists spectroscopic follow up. 
We combine catalogues with spectroscopically confirmed 
LBGs and photometrically selected LBGs which were 
identified as interlopers from 3 < 2 < 6 using observations 
reported by [Vanzella et al.| ([^09[); |Reddy et al.| ( |2006[ ); 
Malhotra et al. (20051; Steidel et al. (20031. We compare 


the fraction of interlopers for LBGs within 5. 0 of bright, 
red foreground galaxies (nir < —22 mag) with the fraction 
of interlopers in the total sample. In both cases, we find 
the interloper fraction to be ~ 8 %, with 2 of 25 LBGs 
around bright foreground galaxies identified as interlopers, 
and 21 of 252 of the entire sample. This indicates that the 
alignment between bright LBGs and massive foreground 
galaxies is not likely due to selection bias. The large 
enhancement in false identifications required to mimic the 
observed magnification bias of ~ 10 for bright galaxies is 
clearly inconsistent with the spectroscopic data. 


7 MAGNIFICATION BIAS AND THE 
LUMINOSITY FUNCTION 

The effect of magnification bias on determining the luminos¬ 
ity function is an important consideration when making a 


census of galaxies in the epoch of reionization {Mason et al. 
|2015[ [Wyithe et al.|2011[ ). In this section we show the effect 
that the lensed fraction reported in this paper has on the 
observed luminosity function. 

The observed luminosity function of LBGs is the convo¬ 
lution between the intrinsic luminosity function, "^{L), and 
the magnification distribution of an SIS, weighted by 
the strong lensing optical depth, r. This results in an ob¬ 
served LF, 4tobs(L), with a power law tail at the bright end 
with a slope of —3 (the slope of the magnification distribu¬ 
tion of an SIS). We assess the affect of gravitational lensing 
on the LF by following the method presented by |Wyithe| 


et al. ( 20111 , where it is modelled by considering the op¬ 
tical depth, r, the mean magnification of multiply-imaged 
sources for an SIS, (/imuit) = 4, and the demagnification 
of unlensed sources (to conserve total flux on the cosmic 
sphere), ^demag = (1 - (Mmuit)T)/(l - t). The observed LF 
is then given by 


'I'obs(L) =(1 - t) -4t(L/^demag) 

/^demag 


+ T 


/■ 


^ dPm,2 




(13) 


where 


d/i 


— 2/(/.i + 1)^ for 0 < /i < 00 is the second 

d-Pm 1 • 


dfi 


image’s magnification probability distribution, and 
given by equation We use the values of the optical depth 
from our analysis presented in Section and calculate 
the optical depth at 2 = 6 . 8 , 2 = 7.9 and 2 = 10.4 to be 
r = 0.72%, 0.80% and 0.94%, respectively. 

We begin by assuming that the observed luminos¬ 
ity function is not affected by gravitational lensing and 
hence represent the intrinsic luminosity functions. We 


plot these intrinsic luminosity functions (Bouwens et al. 
20141, the inferred observed luminosity function, and 


observations (Bowler et al. 2014 Bouwens et al. 20141 at 
2 ~ 4 , 2 ~ 5 , 2 ~ 6 , 2 ~ 7 , 2~8 and 2 ~ 10 in Figure 
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Figure 9. The effect of magnification bias on the bright end 
of the luminosity functions z ^ z ^ 5, z ^ 6, z 7, z 8 
and ^ ~ 10 (LBGs become monotonically more abundant with 
decreasing redshift at Mjjv ^ —22). The observed luminosity 
functions are shown as solid lines, and the intrinsic luminosity 
functions as dashed lines. The luminosity function measurements 
from [Bouwens et al.| | |2014[ | and [Bowler et al.| l |2014| l are plotted 
as circles and diamonds, respectively. At 5: ~ 7, ^ ~ 8 and ^ ~ 10, 
observations are close to probing the bright end where gravita¬ 
tional lensing becomes a significant effect, but not bright enough 
for it to be manifested in the observed LF. 


current surveys out to 2 ~ 8, the affected region of the 
2 ~ 10 LF begins at around Muv ~ —22.5. We investigated 
the 4 unusually bright 2 ~ 10 Ji25-dropouts presented by 


Oesch et al. (20141 to search for evidence of lensing. This 


point is discussed in Oesch et al. (20141, where they find 
there is the possibility of a modest amount of lensing in 
2 of the 4 dropouts. By applying the technique employed 
in this paper, we assign likelihoods of lensing to the four 
2 ~ 9 — 10 LBGs. 

As noted by Oesch et al. (20141, 2 of the LBGs are not 
close in projection to any foreground objects {GN-zlO-3 
and GN-z9-l in the notation of Oesch et al. 20141, while 
the other two do have projected neighbours {GN-zlO-1 and 
GN-zlO-2). GN-zlO-1 is l'.'2 from a foreground galaxy at 
Zphot = 1-6 with Mb = —20.1 (using photometry from the 
3D-HST catalogue). Oesch et al. (2014i infer a photometric 
redshift of 2 = 1.8). Lfsing our redshift-evolving FJR, this 
corresponds to a stellar velocity dispersion of 140 kms“^. 
The required stellar velocity dispersion for strong-lensing 
in this case is 198 kms”^, giving this LBG a likelihood of 
lensing of Jzf = 0.12. GN-zlO-2 is 2'.'9 from a bright galaxy 
at Zspec = 1.02 with Mb = —20.7. This corresponds to 
an inferred stellar velocity dispersion of 214 kms“^. The 
required stellar velocity dispersion for strong lensing is 279 
kms“^, giving a likelihood of lensing of M’ = 0.1. While 
the statistics are too small to draw any firm conclusions, 
this average observed lensed fraction of ~ 6% for the 
four galaxies is consistent with a lensed fraction of LBGs 
brighter than M* of ~ 10%. 


(dashed lines). Figure also shows the biased LFs, 
illustrating the luminosity at which gravitational lensing 
becomes important. This also illustrates the assumption 
that current LF measurements are not significantly affected 
by magnification bias is sound. 

The effect of magnification bias is not significant 
at the faint end of the luminosity function. At around 2 
magnitudes fainter than M*, the excess observed abundance 
of LBGs is on the order of 0.5% for all of the samples, 
which is significantly smaller than the observational errors 
in the abundances at these magnitudes. 

We note that even the brightest Bowler et al. (20141 


and Bouwens et al. (20141 measurements are not bright 
enough to probe the affected region of the LF. However, 
the effect magnification bias will have on surveys at 2 > 8 
is obvious from the solid lines for 2 ~ 8 and 2 ~ 10 where 
the observed LF will display a break from the intrinsic LF 
around Muv ^ —22.5. 

Not plotted in Figure|^are extrapolated LFs at 2 > 10. 
jWyithe et al.j ( |2011| showed that if M* drops sharply at 
high redshifts, surveys of the depth of the XDF with JWST 
will observe galaxies at 2 > 10 in the affected region of the 
LF. 


8 ANALYSIS OF CURRENT J125-DROPOUTS 

In Section [7j we presented the effect that magnification bias 
has on the observed luminosity function. Figure [^highlights 
that while magnification bias is not a significant effect in 


9 SUMMARY 

We have estimated the likelihood of strong gravitational 
lensing of LBGs in the XDF and GOODS at 2 ~ 4, 2 ~ 5, 
2 ~ 6 and 2 ~ 7. We used a calibrated Faber-Jackson rela¬ 
tion to estimate the lensing potential of all foreground ob¬ 
jects in the fields. The result is a measurement of significant 
magnification bias in current high-redshift samples of LBGs. 
Our analysis allows us to draw the following conclusions, 

(i) Approximately 6% of LBGs at 2 ~ 7 brighter than 
M^: (uiHieo ~ 26 mag) are expected to have been strongly 
gravitationally lensed with p > 2; 

(ii) The observed strongly lensed fraction of LBGs at 
all values of mn falls monotonically from 2 ~ 7 to 2 ~ 4, 
which can be explained by the expected evolution in the 
optical depth with redshift, and also M* appearing brighter 
at lower redshift; 

(iii) By evaluating the optical depth in our lensing 
framework, we calculate the magnification bias in each 
sample as a function of M* — Mum, and find that the results 
agree at each redshift and are well described by theoretical 
predictions; 


(iv) Extrapolation of our analysis leads to expectations 
for an increased fraction of strongly lensed galaxies at 
2 > 8, consistent with Wyithe et al. (20111; 
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(v) The magnification bias of the faintest LBGs in the 
sample suggests there may be a flattening of the faint-end 
slope below current detections limits {Muv ^ —16.5). 
However, this result relies on LEG detections at low 
S/N in the XDF, and the constraints are weak. We 
present this result tentatively, with deeper data needed 
to better understand the population of faint high-a galaxies; 


(vi) Assessing the magnification bias as a function of 
luminosity offers an independent method of determining 
Schechter parameters a and M*. The results from this 
method are consistent with those found by fitting the LF 
based on number counts. 


With the confirmation of the role of magnification bias 
come important consequences for future surveys of galaxies 
at 10 < 2 < 20. Gurrently, the LFs at 2 ~ 8 need not be 
corrected for magnification bias. However, magnification 
bias will be significant for luminosity functions at 2 > 10, 
notably in the JWST era ( [Wyithe et al.|[201ip . In par¬ 
ticular, with M* possibly dropping rapidly beyond 2 ~ 8 
(Oesch et al. 20141, JWST will identify predominantly 
gravitationally lensed galaxies at 2 > 10. 
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APPENDIX A: SPATIAL CORRELATIONS 
BETWEEN BRIGHT FOREGROUNDS & LBGS 

In this appendix, we present the manifestation of magnifica¬ 
tion bias in spatial correlations between bright foreground 
objects and bright LBGs, which illustrates the effect 
without relying on the Faber-Jackson relation. 

Source galaxies that have been magnified through grav¬ 
itational leasing are necessarily located in close proximity 
to massive foreground objects. For the leased fractions pre¬ 
sented in Section we expect there to be an excess density 
of bright LBGs around bright foreground objects over the 
average field density. As the lensed fraction decreases with 
decreasing luminosity, the excess probability around bright 
deflectors should also decrease. Similarly, the clustering 
around the more massive, brighter deflectors should be 
stronger than around less massive, fainter deflectors. 

We compute the excess probability of finding an LBG 
brighter than m = 30,28.5,27.5 & 26.5 at z = 7.2 within 
570 of deflectors brighter than some Mb- We choose 570 
as our limit as this is approximately the image-deflector 
separation beyond which strong lensing is unlikely. This 
is confirmed by the distribution of separations shown in 


Figure We also present the spatial correlations for the 
lower redshift samples. At z = 5.9 we examine the same flux 
limits, at z = 4.9 we replace mum = 30 with mum = 25.5, 
and at z = 3.8 we replace miim = 28.5 with mum = 24.5. 
The results are plotted in Figure [Alj 

We find a large enhancement in the probability of find¬ 
ing bright LBGs nearby bright deflectors. As we consider 
fainter LBGs, the excess probability decreases monotoni- 
cally in all of the samples. There exists a considerable excess 
of LBGs around foregrounds with Mb < —23 at ^ ^ 4, 
even for flux limits well beyond M*, however this signal is 
driven mainly by the brightest LBGs. The excess likelihood 
of locating an LBG around a foreground approaches unity 
by Mb ~ —21 for all flux limits in all samples. 

The clustering of bright LBGs nearby massive fore¬ 
ground galaxies is difficult to explain in the absence of 
magnification bias. One mechanism that could produce 
such a signal is the enhancement of LBGs around bright, 
red foregrounds. As discussed in Section |6.3[ we searched 
for this effect in LBG samples with spectroscopic follow-up 
from the literature, and found no evidence that bright, 
red foregrounds enhance LBG detection. Therefore, we 
conclude that the proximity effect shown in Figure |A1| is 
consistent with being due to gravitational magnification of 
background LBGs by massive, bright foreground objects. 
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Figure Al. The excess probability of finding an LBG brighter than various flux limits within 5'/0 of deflectors brighter than Mq at 
2 < 2. In each of the four samples, we And that there is an excess of LBGs around bright foreground objects. The excess becomes 
monotonically more pronounced with brighter flux limits around bright foregrounds in each of the four samples. At each redshift slice, 
we consider a different set of flux limits as M* appears brighter for the lower-redshift samples. The right panel shows a large excess of 
bright LBGs z 7 around bright foreground objects. At 2 ~ 4, we find similar behaviour of bright LBGs appearing more frequently 
around bright foreground objects than in the total field, but the amplitude of the excess is much lower for the same flux limits of LBGs. 
However, for brighter flux limits, we see identical behaviour to that observed at higher redshift. 
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